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Abstract 



We present a new approach for setting initial Cauchy data for multiple black 
hole spacetimes. The method is based upon adopting an initially Kerr-Schild 
form of the metric. In the case of non-spinning holes, the constraint equations 
take a simple hierarchical form which is amenable to direct numerical integra- 
tion. The feasibility of this approach is demonstrated by solving analytically 
the problem of initial data in a perturbed Schwarzschild geometry. 
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I. INTRODUCTION 



The calculation of gravitational radiation from the inspiral of binary black holes is of 
paramount importance to the development of gravitational wave detection into a new form 
of astronomy. The strong field regime of the emission process occurring during the merger 
of two black holes can only be treated by numerical evolution. However, in the near future 
it is unclear whether computing power and techniques will be adequate to track the two 
black holes through many orbits, beginning with a quasi-Newtonian orbit which can be 
approximated by perturbation methods and ending with their merger into a single black 
hole. Instead, the merger radiation may have to be explored beginning with initial data in a 
non-perturbative setting just prior to the final orbital plunge. Numerical codes are presently 
being developed which promise to provide an accurate map between such initial data and 
the emitted gravitational waveform. However, the extent to which the merger waveform will 
have a characteristic signature independent of the details of the initial data is unknown, and 
there is some evidence that the assumptions that have to be made in the usual approach 
lead to data that is not physically realistic [p]]. For this reason it is imperative that many 
different forms of black hole initial data be available. 

Results on multiple black hole initial data were presented by Misner in 1963 ||. As 
is well known he found an analytic solution for the special case of a hypersurface of time- 
symmetry. In this case the extrinsic curvature is zero so the momentum constraints are 
satisfied identically; and, assuming conformal flatness of the initial 3-geometry, the Hamil- 
tonian constraint reduces to Laplaces's equation for the conformal factor. Subsequent work 
has been essentially a generalization of Misner's result. The first step in treating the non- 
time-symmetric case was the development of the theory of conformal decomposition ||. This 
then led to methods for determining initial data for one black hole with specified linear and 
angular momentum [[|. A procedure for "summing" such solutions has been described ||, 
and the procedure has been implemented numerically ||[7|]- Recently, a modified version of 
this procedure has been formulated which avoids complications at the inner boundary for 
this elliptic problem by compactifying the region interior to the black holes ||. However, 
only conformally flat geometries with trace-free extrinsic curvatures can be readily obtained 
by this technique. 

Here we describe the generation of initial Cauchy data by means of a new approach, 
based upon adoption of an instantaneously Kerr-Schild |9|,|10| form of the metric on the initial 
time slice. In this approach the time-symmetric Misner data is not expressible as a simple 
special case. For a single non-spinning vacuum black hole at rest, it leads to Schwarzschild 
initial data in Eddington-Finklestein coordinates. In the numerical treatment of the non- 
vacuum spherically symmetric gravitational collapse of a scalar field, the most stable Cauchy 
evolution algorithms for the inner boundary at the apparent horizon have been based upon 



a generalization of ingoing Eddington-Finklestein coordinates |TTJ. An analogous gauge is 



likely to be beneficial to stability in the two black hole case. Since Eddington-Finklestein 
coordinates are closely related to Kerr-Schild coordinates, the initial data generated by 
our approach offers the additional advantage of allowing greater flexibility in the choice of 
gauge conditions on the lapse and shift in the search for stable numerical algorithms for the 
simulation of black holes. Similar motivations have led to another independent formulation 
of the initial data problem which reduces to the Kerr-Schild form for the Schwarzschild or 
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Kerr cases [112 



The framework used is that of an instantaneously Kerr-Schild space-time, which is de- 
fined and discussed in section O. The corresponding initial data is more general than for 
exact Kerr-Schild space-times, and the constraint equations are derived from first principles 
in section [III]. A numerical procedure for solving the constraint equations is described in 
section |IV|; also, we show that a simple algebraic condition on the metric defines the loca- 
tion of the apparent horizon. We do not, in this paper, solve the numerical problem; but we 
show in section |V] that our solution procedure is feasible by formulating, and then solving 
analytically, the Schwarzschild perturbation problem corresponding to the "close approxi- 
mation" for two black holes |L3] . We conclude in section [VT| with comments on the future 
development of this research. 



II. INSTANTANEOUSLY KERR-SCHILD SPACE-TIMES 



In the 1960s the work of Trautman [H|, and of Kerr and Schild P,|IT|, led to the develop- 
ment of the theory of Kerr-Schild geometries; the theory is comprehensively reviewed in the 
book by D. Kramer et al [ I5fl . Here we use notation such that Greek indices are space-time 



indices and Latin indices are spatial. We will often use a dot to denote a time derivative, 
e.g. V = d t V. 

Let the scalar field V and the null vector field k a define a Kerr-Schild metric 



g a /3 = - 2Vk a kp, (2.1) 

where k a = (—1, fcj) is null. We choose the background to be a standard inertial coordinate 
system with space-time foliated by time slices t = const and the Minkowski metric given by 
Va/3 =diag(-l, 1,1,1). The contravariant null vector satisfies k a = g a ^kp = rj a ^kp. For the 
spatial components we write k a = (1, A:*) so that k l = 5^kj. As an example consider the 
Schwarzschild geometry. In ingoing Kerr-Schild form it is 

M 

V = , kidx 1 = -dr. (2.2 

r 

Another example is the Kerr solution, which is rather long [16[ . 



We wish to pose black hole Cauchy data at t = 0. We only require that initially the 
geometry be of the Kerr-Schild type. In particular, we assume that the evolution of the 
initial data leads to a metric of the form g a p + t 2 j a/ 3, with j a b well behaved at t — and 
9a/3 given by equation (|2.1|). The initial data at t — must satisfy the constraint equations 
G a bn h = 0, where n a is the unit normal to the foliation. These equations do not contain 
second time derivatives so that j a b does not enter and we can analyze them assuming a 
pure Kerr-Schild form, i.e. setting j ab = 0. In a 3 + 1 decomposition, these constraints 
become spatial equations for the initial values of V, ki and their first time derivative, which 
determine the initial values of the 3-geometry and extrinsic curvature. For the family of 
global vacuum Kerr-Schild solutions, the null vector field k a is necessarily geodesic and 
shear-free. The shear-free property makes it unlikely that a globally Kerr-Schild space-time 
can contain more than one black hole. Here, for our instantaneously Kerr-Schild spacetimes, 
we impose no differential conditions on k a except those that follow from the null condition, 
d a (k%) = 0. 
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III. THE CONSTRAINT EQUATIONS 



For completeness, we present the equations in covariant form, and also in "3 + 1" for- 
malism with partial derivatives. 



A. The covariant constraint equations 

Space-time indices take on the values a, j3, 7, 5, . . . = 0,1,2,3. The symmetric part of 
a tensor is T( a/3 ) = ^(T a p + TJg Q ), and the antisymmetric part is T[ a/3 ] = \{T a p — Tp a ). 
The operator for the derivative along the null congruence k a is D = - a k a . The acceleration 
vector of the null congruence k a is a a = k^.^k^ = k a k^ = Dk a . The square of its norm i 



is 



a 2 = a a a a , while the twist of the congruence, u, is defined by u 2 = ^k^.^k ' 13 . 

We may compute the Ricci tensor in Euclidean coordinates explicitly in terms of ordinary 
derivatives as 

Rps = Lp 6 + 2V[(DV)k a a - 2V a a a + D 2 V + 8VW]k p k s 
+4V[2(DV) + Vk a Ja {f3 k s) + 2V 2 a p a s 

+AV 2 k {p k s) ^k a k\ (3.1) 

where 

= v ary (Vkpk s ), aj - {VWkp), 5l - {Vk^k 5 )^ (3.2) 

and D 2 V = V a/3 k a k 13 + V a a a . The Ricci scalar is 

R = -2{Vk a k p ) taP + 2V 2 a 2 . (3.3) 

The Ricci tensor may also be put in covariant form 

Rps = {Vk p k 5 y a . a - (Vk a k s ) ;a p - {Vk a kp), a6 

+2V[D 2 V + (Va a ). a - a 2 V 2 ]kpk 5 + [2V(DV) - V 2 k a . a ]D(kpk s ) 
+V 2 D 2 (kpk s ) + 2V 2 a a k a . w k s) 

and the covariant Ricci scalar is 

R = -2(yk a k^). a(i + 2V 2 a 2 . (3.4) 

Using the expressions for the Ricci tensor and Ricci scalar above we may easily write out 
the Einstein tensor Gp$ = Rps — ^Rgps- The formulae in this section were calculated by 
hand and checked using the computer algebra system MathTensor. 
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B. The constraint equations in "3 + 1" form 



It is straightforward to show that the 3-metric 7^, the lapse a and the shift Pi take the 
following forms: 

lij = ~ ZVkikj, a = —===, Pi = 2Vki. (3.5) 

The ADM constraint equations are found from the extrinsic curvature Kij, which is defined 
by 

Ky = —{Pi\j + Pj\i - jij), (3.6) 

where | denotes covariant differentiation with respect to the metric 7^. The standard form 
of the ADM constraints is 

{3) R + {K\f - K l jK^i = 2a 2 G 00 = 0, K j i}j - (K^)^ = aG°j = 0, (3.7) 

where the first equation is known as the Hamiltonian constraint, and the second equation 
is the momentum constraint; G 00 and G°j are components of the 4-dimensional Einstein 
tensor. In equation ( |3.7| ) indices are raised by the metric 7 U , but we will not be using this 
convention below. 

We have carried out an algebraic computation using REDUCE to evaluate the constraints 
for the instantaneous Kerr-Schild metric. Before we give the results, which are obtained 
after some simplification of the REDUCE output, some general comments are in order. 
The calculation is carried out in the Euclidean background metric where the components 
of covariant and contravariant vectors are identical. Thus it is possible to write all indices 
as suffices, and we choose to do this because it greatly simplifies the computer algebra. A 
repeated latin suffix means summation with respect to the metric b l K These conventions 
will also apply throughout the rest of this section, and in section [IV A|. Some simplification 



rules were used in the computer algebra, explicitly kiki = 1, kikij = and kiki = 0. The 
quantity Zi is defined to be a vector orthogonal to ki, i.e. Ziki = 0; thus Zi is not uniquely 
defined but instead lies in a 2-D surface. The analogue of the 4-acceleration a a is denoted 
by Aj defined by Aj = kjjki. 

The momentum constraint quantity aG°jkj is denoted by , and is 



0= m k = 
-V-2U + 1 



(2V-1) 

-Vijkikj + Vjj + Vk jtj - (4V + l)VjAj + 2A j A j V 2 - /.; ; ; \ '.,/.-, 
-kjdAV - l)V d - 2AjV 2 ) 



(3.8) 



In order to eliminate V, the Hamiltonian constraint is not expressed in the standard ADM 
form. Rather, we use G°o, which is 
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G° = -a 2 G 00 + a2Vm k . 



(3.9) 



Denoting G° by h 2 , we find 

0= h 2 = 
-1 



ViAi + 2kijk jti V - Ak hji kiV 2 - 8A jd V 3 



(2(2V- l) 2 ) 

+4k i j i h i V - 2k hll k 3 V - AkiAiV 2 + \k,. :u k,Y- ) (2V - 1) 

\ki,iY ~\~ 2V,iki)kjj (kijkijV + 1) V^jj 4~ V^jkikj 



(2V-lf 



+2k jd V 2 + kjkjV 2 



+2(2V - l) 2 - ViVjkikj + VjVj - AA^kjhV 2 +Ak lS k l k 3 k l V 2 
+2(41/ -1)(V- tfAjAjV 2 - (32V 3 - 1QV 2 - 8V + 3)VjAj 

-4(W J V-2V j -A j V)(2V-l)k j V 
The momentum constraint quantity aG°jZj is denoted by m z , and is 



0= m. 



Vikj-j - VA t + V, 



^(-2V+l) 
(2V- l) 2 

- fci^-y - ktjkjV ) (2V - 1) 



+2(V j k j ), i V - 2(V r A),i + v ,H k 3 - (hjV + V 1 k j )(2V - l)h 
-[2{V 3 - A 3 V 2 ) - (2V - 1)%F](2V - 



-[2(% + A,)(2y-i)y 2 -y#. 



J:'' 



(3.10) 



(3.11) 



IV. ON THE SOLUTION OF THE CONSTRAINT EQUATIONS 

A. The surface orthogonal case 

In the case of a black hole with spin, such as a Kerr black hole, ki is proportional to 
Vi4> + Q with Ci tangent to the surfaces given by = const. While the analysis of the 
solution of the constraint equations in this general case is quite complicated it simplifies 
drastically if q = so that ki is orthogonal to 2-surfaces S, with 



(4.1) 
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The constraints may then be solved by carrying out a hierarchy of integrals along the rays 
normal to the surfaces S. (Note that difficulties may arise at points where Vj0 = 0). 
Here, we confine our attention to this surface orthogonal case. It includes the description of 
boosted black holes. 

Formally, the ansatz used is that is given everywhere on a 3-D manifold with Euclidean 
metric 5ij. The vector field ki is defined by equation (|4.1|) , and is everywhere orthogonal 
to 2-surfaces 5* (which are equipotentials of <p). Also, V and ki are given on a particular 
equipotential So- The solution procedure assumes that ki is defined everywhere, and also 
that ki t i 7^ everywhere. In multiple black hole problems it is likely that these conditions 
could be violated at isolated points, and special procedures would need to be developed for 
integrating through such points. 

The equations h 2 = 0, = and m z = appear complicated, but it should be 
remembered that we know terms involving k{ and any spatial derivative, e.g. A^. In 
addition, we know terms involving V, ki and any spatial derivative within Sq. In particular, 
the following terms are known on So from the given data: 

• ViAi and Vjkj. Since both Ai and kj are orthogonal to ki they both lie in So, and 
therefore the V^-derivative in both terms is calculable from data on So- 

• Vjj — Vijkikj and —kjVjkiVi + (Vj) 2 . In both cases the off- So derivatives cancel out, 
and so the terms are calculable from data on So- 

• ki t i. Since ki is orthogonal to ki, it follows that ki t i can be expressed in terms of 
derivatives of ki within S . 

Equation (|3.10|) (h 2 = 0) is a linear equation for kiVi, and may be written schematically as 

d k V = 1 {D 1 ), (4.2) 
where O represents an operator on S , and D represents the data given on S ; d k means 



kidi. Next, V is found explicitly from equation (3.8) (m& = 0): 



V = 2 (D 2 , d k V). (4.3) 

Finally, equation (|3.11|) (rn z = 0) is regarded as a linear equation to determine kijkjZi 
(Noting that Zidi means a derivative within So, it is straightforward to see that all terms 
are expressible entirely in terms of known data): 

z t d k k t = 3 (D 3 ,d k V,V). (4.4) 



Because Z{ can take two independent directions, equation (fOJ) gives two components of 



remaining component follows from the condition that fcjfcj = so that k^jki 



jki, cincl ttiGrcforG jkjk^ — k^ jkjki 



In this way V and the d k derivatives of V and kj are found at each point of So. Thus 
the solution may be extended to the "next" 2-surface S, and, in the absence of singularities, 
to the whole manifold. 
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B. Multiple black holes 



For a single (unboosted) Schwarzschild black hole, fcj is a radial field which can be ob- 
tained by the equipotentials of the function = 1/r. Thus the choice = EjMj/rj generates 
Schwarzschild-like equipotentials about each site where r\ = and becomes a candidate for 
data for (unboosted) multiple black holes with masses Mj. By giving appropriate data for 
V on So we shall show that this does produce the necessary apparent horizons. 

In order to describe boosted black holes we begin with the null vector k a = —V a (t + r) 
for a single Schwarzschild black hole. Under a boost with velocity v in the z-direction, 
t — > (t + vz) I y/l — v 2 and z — > (z + vt)/ y/\ — v 2 , we have 

t + vz+y/(l- v 2 )(x 2 + y 2 ) + (z + vt) 2 

t + r " VT^ • (45) 

Then, for a single boosted Schwarzschild black hole, initial data for fcj at t — can be 
obtained from the equipotentials of = {vz + [(1 — v 2 )(x 2 + y 2 ) + z 2 ] 1 ^ 2 }^ 1 

Construction of a solution corresponding to multiple black holes requires a numerical 
integration procedure. For simplicity, we restrict our discussion to an analysis of the binary 
case <fi = Mi/ri + M 2 /r 2 , although the above considerations make it apparent how to 
generalize to the case of multiple, boosted black holes. In this simple case, the equipotentials 
consist of small disjoint (topological) spheres which merge in a "figure-eight" to form a 
common set of outer spheres. The hard problems arise at the bifurcation point of the 
"figure-eight". We could avoid this problem by considering the case where the two holes 
are close so that they are surrounded by a common apparent horizon. This would suffice 
to determine data for the exterior spacetime without going down to the bifurcation point. 
The problem with doing this is that we would probably not be able to include the merger 
radiation this way and that is the most interesting feature. Nevertheless, it seems advisable 
to tackle the close binary case first. 



C. Apparent horizons 

In order to make sure that the data proposed above really represents black hole data we 
must demonstrate the existence of apparent horizons. Consider the equipotential 2-surfaces 
S given by = const. For the reasons given in the introduction, we are interested in ingoing 
Kerr-Schild coordinates so that k a is chosen to be the incoming null vector normal to S. The 
divergence of k a is the same when measured either in the full Kerr-Schild metric or in the 
background Minkowski metric. Thus it is apparent from the asymptotically spherical shape 
of the equipotentials on the approach to each r\ = site that p- m (the incoming divergence 
of S) goes uniformly negative. 

The outgoing null normal to S is given by £ a = 2t a — (1 + 2V)k a , where t a — (1, 0, 0, 0) 
is the time translation Killing vector of the Minkowski background. It has normalization 
£ a k a = —2. Let (£ a , k a , m a , rh a ) form a null tetrad. We choose normalization m a m a = 2 
(with respect to both the Kerr-Schild and Minkowski metrics). The complex spatial vector 
m a is tangent to S. The outward divergence of S is given by p ou ^ = m' a m^V a ^. We have 
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2m (Q m /3) V Cf t /3 = m {a Th p) C t g a p = m^fh® X t r] a p = 0, 



(4.6) 



since t a is a Minkowski Killing vector. As a result, p ou t = —(1 + 2V)P[ n - 

Where |V| is small, p- m is negative and p ou ^ is positive, as for a sphere in Minkowski 
space. For the Schwarzschild solution V = —M/r and the (apparent) horizon forms at 
V = —1/2. Our result shows this holds in general: An equipotential with V = —1/2 is 
marginally trapped and determines the position of the apparent horizon on the initial time 
slice. 

Thus, since V is free data on some choice of equipotential Sq (which might consist of 
disjoint pieces), it is straightforward to give data that guarantees the existence of black 
holes. Note that the location of apparent horizons is considerably simpler here than in the 
standard approach. This suggests that the apparent horizon itself be chosen as the inner 
equipotential to start the integration process. 



V. PERTURB ATIVE CALCULATION 

We use the above methods to find initial data for a surface orthogonal problem, specifi- 
cally a perturbed Schwarzschild geometry. 



A. The close approximation 

As described in Sec. [IV B| , two non-spinning black holes initially at rest can be modeled 
in terms of the potential = M\/r\ + M 2 /r 2 , where T\ = \x l — x\\ and r 2 = \x l — x\\. 
Without loss of generality, we can take the line connecting the black holes to be the z-axis, 
with the origin located according to the center of mass condition M\Z\ + M 2 z 2 = 0. We let 
a = Z2 — z\ and M = M\ + M 2 , and use standard spherical polar coordinates (r, 9, (p). The 
close approximation is defined by the condition 

« = s « I- ("I 



We find the expansion 



with 0(e 3 ) remainder, where 



M e 2 M 1 M 2 MP 2 . . 

+ (5.2) 



P 2 = (3cos 2 fl-l)/2. (5.3) 



The approximation is valid when r > max(zi, z 2 ), and improves with increasing r. To this 
order of approximation, the covariant components of k in standard spherical coordinates are 

ft. t n ( 1 3e 2 M!M 2 costing ^ 

[k r , k e , k<p) = (-1, , 0). (5.4) 

The ansatz on the initial hypersurface t = is then equation (|5.4p above, together with 
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V = -M/r + e 2 v(r)P 2 



(5.5) 



V = e 2 v T (r)P 2 



(5.6) 



(k r ,kg,k<p) = (0,3e 2 k T (r) cos 9 sin 9, 0). 



(5.7) 



The functions vt{t) and £/r( r ) are, in a 4- dimensional sense, time derivatives; but since the 
whole analysis is carried out in a spacelike 3-D manifold, formally we regard Vt(t) and kr{r) 
as independent functions on the manifold rather than as derivatives. The objective now is 
to determine the functions v (r), Vt(t) and kx{r). 

B. The algebraic computations 



The computation of the constraint equations has been done in two different ways, and 
using two different algebra systems (Maple and REDUCE). In one calculation we put the 
ansatz above into the equations reported in section [III B| , and in the other calculation the 
ansatz is put into the full 4-dimensional Einstein equations and G° a is found. Both pro- 
cedures give identical results, and this serves as a partial confirmation that our computer 
algebra programs are correct. 

The constraint equations = 0, h 2 = and m z = were evaluated, and to order zero 
in e they were satisfied identically; and to order e 2 the angular terms cancelled out and the 
equations involved r only. After some manipulation we find: 



6k T (r)M 2 + 3Mrk T {r) + 



6MiM 2 M 



{r - M} + 3v (r)r 2 



(5.8) 



v'{r) 



Av (r)r 2 + 



§M X M 2 M 



{r-M} + 6k T (r)M 2 



(5.9) 



k' T (r) 



Mr 2 



2Mrk T (r) 



6MiM 2 M 



v(r)r 2 



(5.10) 



where ' represents Equations ( |5.9|) and ( |5.10|) constitute a system of two first-order equa- 
tions for v'(r) and k' T (r); once they have been solved vt{t) is found explicitly from equation 
(|5.8|). The system ( |5T9l) and ( |5.10| ) may be re-expressed as one second-order equation: 



v"(r)r 5 + 5rV(r) - Qr 2 Mv(r) - 6M i M 2 M | 3r + 2M } = 0. 



(5.11) 
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C. Behavior of the solution 



In this section we are concerned with finding the behavior of the solution of equation 
( 5.11]) . Making the transformation 



24M 



r = — — , v(r) = x u(x), (5-12) 
x 2 



equation (|5.11| ) becomes 



x 2 u, xx + xu iX - u(x 2 + 16) - 6 ^2^2 ( 36a;2 + x4 ) = °- ( 5 - 13 ) 



The solution to equation ([5 .13 ) is 



/ \ t / \ mr/- \ M\Ml, , , X 2 lJx)\0gX, i . 

u{x) = CJ^x) + C 2 Y A (tx) + -j^(w(x) - 2304 - ^ ), (5-14) 



where I4 and I4 are standard Bessel functions (T^| and i = v^T- The function w(x) may 
be written as 

w{x) = wqx 6 + wgx 8 + wiqx 10 + ... (5.15) 

and the coefficients in the series are found from the following recurrence relation applied to 
n =6, 8, 10..., and using the initial condition iy 4 = 0: 

[n 2 - 16K = w n . 2 + 36( n_ 2)!(i + 2)!2w - (5-16) 



Using equation ( |5.12| ) we may now find v(r); then equations ( |5.9|) and ( p.8| ) lead to k T (r) 
and Vt(t). 



The leading term in the inhomogeneous part of the solution (|5.14|) behaves as x 2 as x — > 0, 



which means that v(r) would behave as r -3 as r — ► 00; this is the expected quadrupole 
behavior at infinity. The constants C\ and C 2 represent gauge freedom, as is evident from 
the fact that they persist in the Schwarzschild case (Mi = M and M 2 = 0). I4 behaves as 
x~ 4 near x = 0, so that if C 2 7^ then v(r) would behave as a constant as r — ► 00. While 
this is possible, equation ( |5.9| ) then indicates that kxir) would be singular as r — > 00. Thus 
it is convenient to fix C 2 = 0. I4 behaves as x 4 near x = 0, so that v{r) would behave as 
r~ 4 as r — > 00. C\ could be chosen so as to fix the position of the apparent horizon: for 
example, C\ could be chosen so that v (2M) = 0, which would mean that the position of the 
apparent horizon is not perturbed relative to the Schwarzschild case. 



VI. CONCLUSION 

In this paper we have shown how to express the initial value problem of general relativity 
in a way that can be solved using initial data that is instantaneously Kerr-Schild. We 
have used this approach to formulate, and solve analytically, the initial value problem of a 
perturbed Schwarzschild geometry. It may be interesting to see if the difference between 
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this initial data and that obtained from the conformally flat method affects results about 
the perturbative evolution of a black hole [13]. We have shown that the apparent horizon is 



described in a very simple way, V = — |; thus the generation of initial data for black holes 
by the Kerr-Schild approach would automatically give the initial location of the apparent 
horizon. 

Clearly the Kerr-Schild approach to the initial value problem is very much in an early 
stage of development, particularly in comparison with the conformally flat method that was 
introduced in 1963 0. This paper is just a first step towards developing this new approach. 
The next steps, which are currently under investigation, will be the construction of a code to 
solve numerically equations (|3.8|), fl3.10| ) and ( |3.11| ), and the generalization of our approach 



to be able to include problems with spin, where hi is not orthogonal to the 2-surfaces S. 
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